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After presenting a brief review of the phenomenology of the Leading Effect, we define a new 
variable, the "Total Available Quadri-Scalar" (Atot) and propose it as the invariant quantity 
effectively available for the production of the multihadronic final states. The introduction and 
V) ' definition of this new varaible are justified by means of simple geometrical-kinematical consid- 

erations and we show that ^4tot reduces to the so-called effective energy in the single specific 
situation where the use of the latter applies. Using Atot to re-plot existing data, the quantity 
(ri c h) is shown to be a "Universality Feature" - that is, independent from the colliding particles, 
the collider nominal energy, and even from the hadronic invariant mass - as imposed by QCD 
universality. 

a: 

^ ■ 1 Introduction 



Basically speaking, QCD shows two main features. One of them is perturbative in nature and 
is Asymptotic Freedom. This property is understood in terms of a negative j3 function, which 
implies that the QCD gauge coupling as decrease with an increasing q 2 . The other, Confmament, 
is non-perturbative and has not received, up to now, a satisfactory theoretical explanation. 

Another important but almost forgotten feature of QCD is the Effective Energy. This is also a 
non-perturbative effect and is, roughly speaking, the mechanism through which the initial (nom- 
inal) energy is shared among different processes, one of which is the hadronization. The slice of 
energy that goes to the latter is not the nominal, but an effective energy. It is this one, and not 
the whole nominal energy, that is at disposal for particle production, as explained in the paper in 
a more datailed way. 

This very interesting property is practically no more studied, even if it could shed some light on 
and provide an alternative approach to the question of how the hadronization mechanism works. 
The fundamental idea that lies behind the effective energy approach is that of distinguishing two 
main phases inside any given interaction, namely, the quantum number conservation (or flow), 
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and the hadronization in the proper sense of the word, that is the process by which quarks and 
gluons "hadronize" and become observable matter (hadrons). 

Some papers published in the early '80s [1, 2, 3, 4] showed that, whatever the invariant quantity 
available to the hadronization process, it should equate the total hadronic energy of the whole 
system as evaluated in the CMS (therefrom the name effective energy). 

In fact, if a whole set of quantities (specified in the next section) are studied in terms of 
that portion of energy, the plots relative to different kinds of collisions between different kinds of 
particles converge toward the same curve. 

Despite the great importance of such a result, neither the correct invariant representation has 
been identified, nor a satisfactory justification for the choice of the quantities adopted in the past 
to represent it has been given. 

In fact, the quantity commonly used to plot data and to study the world- widely collected 
experimental results is Mx = y (<Ztot) 2 but, as will be showed in this work, this quantity cannot 
work if we insist in separating the aforesaid phases of the interaction. Furthermore, as some 
recent papers show, (see e.g. [5]) M x is unable to reveal the Universality Features when also DIS 
processes are taken into consideration. 

Purposes of this paper are to introduce, by means of simple kinematical considerations, what 
is believed to be the correct Lorentz-invariant representation of the quantity from now on called 
the "Available Quadri-Scalar" and to show how it is possible, by correcting the above mentioned 
result, to let the universality features be "revealed again". 

The outline of this paper is as follows: in Section 2 a brief review of the phenomenology of 
the Leading Effect is presented; in Section 3 a new variable, ^4tot , is conceptually introduced 
while its formal definition is given in Section 4. In Section 5 the relations between Atot and 
Mx are examined. Evidence of (n c h) vs vItot universality is provided in Section 6 with the final 
result given in Section 7. In Section 8 rescaling of p-p data is discussed and in Section 9 a brief 
analysis of the dependence of our result from kinematical condition is given. Section 10 contains 
the conclusions. 

2 Phenomenology of the Leading Effect 

The fact that the total energy available to particles production in a given type of interaction is 
not, in general, the nominal energy of the reaction, but another quantity that takes into account 
the leading effect, is one of the most important discoveries made in the early '80s [6]. Before 
that time, all the measured quantities were analysed in terms of y/s and this brought to different 
results in different experiments. This was commonly accepted even if it was in a flagrant contrast 
with the QCD universality. 

In fact, at a fundamental level, the final state of any interaction should depend only on some 
(Lorentz-invariant) scalar variable believed to be available to the hadronization mechanism, but 
not on how that quantity has been put together. This means that, had this quantity been yfs , the 
results of any analysis should have been the same, independently on the kind of reaction under 
exam (e.g. e + e _ , p-p , DIS). Instead, the results were all different and there was no explanation 
for this situation, to which people referred as "the hidden side of QCD". 

As mentioned, in the early '80s it was pointed out [1] that QCD universality could have been 
made explicit if the quantity E™ , that is, the total hadronic energy of the reaction given by 
subtracting the energy of the leading particle(s) from the nominal energy of the reaction, was 
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used to plot data instead of \fs . Within a few years, this effect, called the Leading Effect, 
would have been shown to be universal: in fact, no matter if the interaction studied was strong, 
electromagnetic or weak, the leading effect was always present [6, 7, 8, 9]. 

was soon after called the "Effective Energy", the name claiming for -E"ot to be the portion 
of energy effectively available to the production of the multihadronic systems present in the final 
state, after subtraction of the energy of the leading particle. 

The leading particle was defined to be the particle leaving the interaction vertex with the 
highest longitudinal momentum [1]. The role of this particle was to carry, totally or partially, the 
quantum numbers (as J PC or flavour) from the initial to the final state. The transfer of these 
quantum numbers from the reacting particle (s) to the leading particle was called the Quantum 
Number Flow (hereafter: QNF). 

From 1980 to 1984 a series of experiments were conducted at the ISR (CERN) to establish if 
£™ was the effective energy. All these experiments proved that there were no differences between 
p-p and e+e~ collisions results if E"*° was used to perform the analyses. The quantities measured 
in these experiments were called the "Universality Features" as they showed the same behaviour 
whatever the kind of experiment and the nominal energy of the collider. Some of them are: 



1 ) (n c h) — Mean charged particle multiplicity [3] 

2 ) = Fractional energy distribution [1, 2] 

3 ) -—J = Transverse squared momentum distribution [10] 



dp: 



4 ) — — ^ — r = Reduced transverse momentum distribution [11] 

d (<p t >> 

E c h 

5 ) a = ° = "Charged" energy [12] 

-^HAD 

6 ) , „ , = Event planarity [4] 

d(pPp) " ' ' 

\ l 1 I IN, OUT 

7 ) P(n c h) = Charged particle multiplicity distribution [13] 

8 ) iV° of propagating quarks vs L [8] 



Some studies regarding DIS processes were also made, but only using DIS variables to plot p-p 
data, being impossible to re-analyse DIS data in terms of E"£° [14, 15]. 

They showed no more that some agreement between the two curves but were almost useless to 
decide whether E"*° worked in DIS as well as it did with other processes, as the variable used 
was W 2 that does not take into account the leading effect. 

In 1984 the ISR closed and no further intensive experimental work was planned in this important 
field of non-perturbative QCD. The last relevant result was obtained in 1997 at LEP (CERN), 
when the meson r/, very rich in gluon content, was seen to be produced in gluon induced jets as a 
leading particle, that is, having an anomalous high longitudinal momentum [16]. This completed 
the series of experiments aimed to establish the universality of the leading effect in p-p and e + e~ 
processes but nothing of conclusive had been obtained in relation with DIS processes. 
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As invariant representation of £™ was chosen Mx = \J ((7™?) 2 [6> 13, 14, 17] in consistency 
with the previous use of -y/s : again a total invariant mass had been chosen to be the fundamental 
quantity from which the hadronization should have been depending. Indeed Mx reduces to -E"ot 
when evaluated in the CMS and, as all the experiments were made in balanced colliders (where 
CMS=lab), then )cms = (E™ )lab held, so the latter was in effect the right quantity to 

be measured. 

Nowadays, because of an unjustified extention of the use of Mx sto DIS processes, the almost 
totality of the experimental works use Mx to study the behaviour of a given quantity, and the 
invariant mass is generally but wrongly believed to be the correct representation for the effective 
energy. 

As a result, papers about this topic often disagree when they try to establish if, for instance, 
(n c h) has a universal behaviour. See, as an example, [18, 19], where the charged particle multi- 
plicity is analysed. They both agree with the hypotesis of the universality features, also referred 
to as "fragmentation universality". 

On the contrary, in a more recent work [5], a discrepancy of 15% is observed between DIS and 
e+e~ , p-p mean charged particle multiplicity. 

Strangely enough, the use of Mx is going on even when a fundamental quantity as (n c h) has 
been shown to be no more "universal" if also DIS data are considered. 

From the viewpoint adopted here, in [5] is proved that Mx cannot represent the invariant 
quantity effectively available for particle production and here is a first indication that should have 
long since been considered: the logical step that brought to the introduction of the effective energy 
was to separate the interaction in the two already mentioned sub-interactions. Now, if we insist 
on associating the total invariant mass to any given particle system, we should also make the 
association: 

Leading Effect < — ► M L = ^/(q^ D ) 2 

But, as it is easily seen, an invariant-mass-type quantity cannot be considered as the variable to 
be associated to any sub-system into which the whole system is being divided as 

sfs £ M x + M L (1) 

So we have to make our choice: either we separate the final state particles into leading and 
hadronic, or we use invariant masses, but not both. 

3 Introduction of A TO t 

Firstly, let us recall that in the "leading approach", the interaction is divided into two processes: 

Leading Effect = the quantum number conservation mechanism 

Hadronization = the mechanism that transforms some available quantity 
into particles masses 

The quantum numbers of the incident particles flow from the initial to the final state thanks 
to the leading particles (usually but not always identifiable with the so-called remnant) [8] and 
actually carried by them. Because of this "enhanced" dependence from the incident particles, the 
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leading particles usually acquire an anomalous large longitudinal momentum: this is the Leading 
Effect and these particles are said to be produced as leading. 

Then comes the hadronization, that converts what is left into the whole of the other particles 
produced in the interaction. The leading particles and their 4-momenta are no more considered 
in the sense that they are not viewed as part of the hadronic system produced. What is left after 
this subtraction is what is eventually studied and analysed. This is the "effective energy approach" 
presented in few words. 

And here is the issue: is it after all correct to ignore the whole leading particle as it has been 
done so far in all the works published in this field? The answer is: absolutely not. In fact the 
leading particle leaves the vertex following a trajectory that is not perfectly aligned to that of the 
incident particle. For example, if the leading particle is the proton remnant, it is indeed true that 
it will have a large longitudinal momentum, that is, it will be strongly aligned with the incident 
particle axis, but not exactly nor completely. 

Keeping in mind that there are some exceptions to the situation we are about describe, let 
us try to show the problem by schematizing the incident and the leading particles in a generic 
collision, see Figure 1. 



in: 



P 2 



-> LEAD -> LEAD 

Pi P2 



out: 




Fig. 1 Simplified sketch of a collision: in the initial state (top) are represented the two incident particles, 
while in the final state (bottom) only the two leading particles are shown. The transverse contribution to 
their 4-momenta cannot originate from the 4-momenta of the incident quarks and must be a consequence 
of the collision, that is, of the quark(s) loss/exchange: for this reason it cannot be associated to QNF and 
must not be subtracted to y^s (see text). 



In the "out" sketch it is evident that the transverse component of the leading momenta cannot 
originate from the incident momenta in a direct way. In other words, the transverse leading 
momenta are not those of the spectator partons (as implicitly assumed in previous works on the 
leading effect) but represent the partons exchanged in the interaction. 

This remark is three-dimensional, but it is straightforward to introduce the appropriate four- 
dimensional version by simply extending the formalism, that is, decomponing the energy into its 
longitudinal and transverse component. 

As our starting point has been the QNF, and as these quantum numbers ought to be carried 
by the leading particle, it is correct to refer to the "motion" quantum numbers by visualizing the 
trajectory of the leading particles. 

Now, if the longitudinal 4-momentum of the leading particle can be safely associated with that of 
the quantum numbers propagating from the initial state, this cannot be the same for the transverse 
component. The latter must come from partons exchange with the other incident particle, when 
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the leading particle is the same as the incident particle, or from the partons loss/exchange, in the 
case where the leading particle is different from the incident one. 

This simple "geometrical" remark immediately suggests that: 

The component of the gj;™ that is transverse to the QNF is not leading. As such, it 
must not be subtracted from yfs to get the Total available quadri-scalar. 3 

The longitudinal component of will be called £tot and represents the "cost" to be paid 
in a given event to force the quantum numbers conservation. What is left from y/s is the invariant 
variable to be associated with the multihadronic final state and that represents the "slice" of yfs 
effectively available to its production. In other words, it is the Lorentz-invariant quantity that 
should be used to analyse and describe every experimental result. Summarizing we have: 



The indices "tot" stays for "total" and reminds us that the whole product of the interaction 
is being considered. Anyway, in principle, nothing would prevent us from considering one hemi- 
sphere only by defining A\ for hemisphere 1, and A% for hemisphere 2. It is easy, in such case, 
to change the definition of Atot given below: it suffices to consider the leading (or hadronic) 
quadrimomentum measured in one hemisphere instead that looking at both (see next section for 
a consistent definition given in equation 7). This means that it is possible to use the variables 
introduced later in the paper to study target experiments too. Useless to say that this agrees, and 
is imposed by, the universality of the leading effect. This situation will not be analysed here and 
when we talk about the available quadri-scalar we are actually talking about the total available 
quadri-scalar. 

It is important at this point to distinguish two kinds of leading effect that will be called here 
direct and indirect. The first type is observed in all the collisions but the annihilation processes. 
The second one regards the annihilation processes only. 

But how is it possible to have a leading effect if the incident particles annihilate? Actually, 
it is still possible to talk about a leading effect but it relates no more to the propagation of the 
quantum numbers of the incident particles. Rather, it is referred to the quantum numbers of the 
particles formed by pair production from the j/Z° produced at the annihilation vertex. 

This second type of leading effect has been observed in e + e~ annihilations (it is present at a 
1% rate) and the most important studies regarded the production of the charmed meson D* [6] 
and that of the "gluonic" meson rf [16]. 

As for the D* , the propagating quantum numbers are those of a c or c quark while, when the 
rf is concerned, its strong gluonic component tells us that it is carrying the quantum numbers of 
the gluon: in fact, the 7/ appears to be leading only when emerging from a gluon induced jet (3 
or 4-jet events). 

Why this distinction? Apart from the conceptual, there is a practical reason to insist on it. 
In fact, the calculations performed later cannot be referred to the indirect leading effect, and the 

3 Here a short note about the name chosen for this new invariant is appropriate: it is simply the most compact 
among the correct names. In fact, it would be wrong to use the word "energy", as well as it is not correct to 
use the term "mass". In fact Atot is not an invariant-mass-like quantity. So, the good practice to call things 
with a name that resembles their nature has been applied. 




£tot § oes to the leading effect 



-<4tot goes to the hadronization 
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given definition of >1tot itself ought to be changed. It would not be a difficult task, and the 
appropriate formulation in this case is also given below, but the experimental evidence has not 
been pursued yet. 

It is evident anyway that, even in such situation, some projection must be done. We can guess 
that the axis on which q hEAD should be projected is the axis given by the mean momentum of the 
jet containing the leading particle, namely, the "charmed" jet or the "gluonic" jet respectively. 

Keeping in mind this distinction, it is now possible to find a mathematical expression for the 
scalar quantity available to the energy-into-mass transformation. 



4 Formal definition of £ T ot and A TOT 

Following the remarks made in the previous section and using the minkowskyan metric instead of 
the euclidean, it is immediate to write the expression for £tot 

INC 

r „LEAD 'TOT (c\\ 

i-TOT = <7 TOT ' —T \ l ) 

V s 

According to our point of view, this represents the cost of the quantum numbers conservation 
in a given event. What is left is what the physics assign to particle production 

^TOT = V(<^-C A T D -^f (3) 

or, in compact notation 

= V* — -Ctot (4) 

This is the correct formal definition of A TO t but we should now find a manageable expression to 
be used below. A possible one is given by writing 

INC \2 



V^F=^§= (5) 
and substituting in (3) to get 

INC 

A HAD ''TOT /£\ 

^TOT — 9tot ' —r \°) 

V s 

that is, ^4tot is the projection of the sum of all the 4-momenta of the produced hadrons (excluding 
of course the leading particles) on the axis given by the incident particles. 

Incidentally, this last remark suggests the only possible consistent definition for A\, that, as 
discussed in the previous section, is the variable to be used when considering a single hemisphere 
(in which case also A2 could be of interest) or a target experiment: 



INC 

A l = q^ D ■ (7) 
V s 

It must be recalled again that this result is only valid for the direct leading effect. When we 
are concerned with the indirect leading effect we should write 

£ TOT (ind) ee (8) 
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and coherently modify the expression for ^4tot • The last is valid provided that there is no further 
leading effect observed in the other jets (in which case the definition of £tot should be changed 
accordingly) . 

First of all, the consistency with previous works must be tested. These all clearly showed that, 
if the variable 

R HAD R INC _ PLEAD /q\ 

^TOT ^TOT ^TOT 

is used to perform data analysis, then universality features manifest themselves: a whole set of 
measured quantities shows the same behaviour independently of everything but -E"ot • So the 
value of A TO t in terms of E"£° must be estimated. 

The colliders where the leading effect was studied (ISR, LEP, HERA(e + e - )) were all balanced, 
that is, the incident particles had the same energy. In this case 

^tot = (£ t H ot)lab = (Qcms (balanced colliders) (10) 

must hold. In fact, what we get if we specialize to the CMS is: 



~inc / 771INC . n\ 

A /7 HAD • TQT / TTiHAD . — *HAD \ l^TOT' U / c 

^TOT Htot ' r~ v tot > I 'tot/ c 



f— V TOT ' r TOT/ CMb f~Z 

V s V s 

— \-Cj tot )CMS 



(11) 



TOT ' TOT 

/?INC 

TOT / CMS 



and there are no consistency troubles as we recover the variable used in that early works performed 
in the CMS limit. The previous equation also shows that: 

NON CMS 

(12) 

that is a very important result and tells us that, in general, E™ is not the effective energy. In 
other words, at an unbalanced collider, the direct use of -E^ot w ^ n °t work and is not correct. 

This means that the universality features were discovered thanks to the fact that the exper- 
iments were incidentally conducted in the system where the available quadri-scalar matches the 
energy of the hadrons produced (excluding the leading particles). Probably, this is also the reason 
that induced the discoverers to think that a non-invariant quantity like the hadronic energy, could 
have had such a fundamental role. Had those colliders been unbalanced (as is HERA nowadays), 
the discovery of the universality features could not have been made. 

In order to be exhaustive, it should be said that some late works made at the ISR were performed 
using Mx instead of -E^ot [6, 13, 14, 17] but the Universality Features were evident anyway, even 
if, in general, Mx ^ E"£° . How is it possible? In chapter 8, it will be showed how the reason 
for Mx worked so well in those occasions, lies in the particular cuts made on the data set at that 
time as well as in the low energies and resolutions. 
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5 Relation between Mx and A TO t 



As mentioned, in general, ^4tot 7^ Mx ■ The best way to get a useful relation between these two 
variables is to estimate the quantity (^xot — ^x) : 



a 

A 2 — C,/e — „ LEAD . TOT x2 _ 9 LEAD INC , 

^TOT — IV s itot 1-) — s z 9tot Vtot + 



/ -yLEAD /ylNC *\2 



/If 2 _ / INC _ LEAD \ 2 _ _ 0„INC LEAD , / LEAD \2 
X V iTOT ^TOT / "iTOT ^TOT ' V^TOT / 

so that the difference is all in the last terms: 

/-LEAD^INC \2 

Wtot "tot/ 



(^TOT ~~ ^X.) 



/ LEADN2 
V^TOT / 



(13) 
(14) 

(15) 



This is the correct invariant expression for the difference we are estimating, but we also need 
an expression that is function of some measured quantities in order to perform the calculations 
that follow. 

To simplify the notation we introduce the following shortcuts for the indices Total, Incident, 
Hadronic and Leading: 



INC 
LEAD 
HAD 
TOT 



I 
L 
II 

T 



Calculating in the CMS we get 



A 2 



Mi 



L\2 

CMS 



- (4? 



L\2 

CMS 



- {P L T ) 



L\2 

CMS 



L\2 l 

CMsJ 



(4? = 

= (ri)Ls = (P 'r ' 



H\2 

CMS 



(16) 



This also follows by observing that 



A 2 

M 2 = 



~ {Et )cms 
( E t)Ims ~ (PT )cms 



(17) 



so that 



Atot ~ M x — (Pr)ci 



(16') 



The importance of equation (16') is threefold: 

1. It shows that ^4tot > Mx that is, the hadronic system "has more 4-scalar at disposal" for 
particle production than what was believed so far. This implies that more particles can be 
produced than the value of Mx could suggest (and at constant Mx) , depending on the total 
3-momentum of the leading particle(s) (see point 3 below) as evaluated in the CMS; 
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2. It gives a very simple recipe to re-analyse data that were wrongly plotted vs Mx: it is in fact 
sufficient to substitute (event by event) Mx with the sum in quadrature of Mx and (p?) cms'- 

Mx — ^ x + (p1)L or (18) 

^tot = ^ m x + (Pt)?ms (18') 

Here the knowledge of (Px)cms ls necessary, otherwise the previous formula becomes useless 
and another method, to be showed later in the paper, must be employed; 

3. It highlights the importance of the unbalance in (p-t)cms between the two hemispheres: what 
makes the difference between ^tot and Mx has nothing to do with the absolute value of the 
leading effect but depends only on its unbalance: the more this unbalance is, the more Mx 
differs from Atot, the less Mx is able to highlight the universality features. In other words, 
you can have a strong, a medium, or a small leading effect: if the momenta of the two leading 
particles detected are the same, then you will always find that Atot = Mx , regardless of 
their absolute value. 

A TOT =Mx = (^) CMS (if (pf +pf) CMS = 0) (19) 

In view of the application of ^tot to compare p-p and e + e~ data, another fact must be 
discussed here, not only for future reference, but also to check the compatibility with previous 
results. 

When the confrontation between p-p and e + e~ data were made, the latter were plotted in 
terms of \fs . When we do not have a leading effect in e + e~ annihilation, the relation y/s =Mx 
holds. But, as stated before, a strong leading effect in e + e~ processes is only present at a 1% 
rate. In other events the leading effect can be neglected. 

Besides, when an average among all the events is performed, the 1% showing a large leading 
effect is easily overcome by the remaining 99% and becomes undetectable, in the sense that the 
difference between Mx and y/s is far below the experimental uncertainties. This means that, for 
those samples of e + e~ data, you have 

sfs = Mx 

and, as e + e~ colliders always have the feature lab=CMS, the relation 

Mx — -Bto? 

holds too. In other words, for our purposes, in e + e~ annihilations: 

A TO t = M x = Vs = (e+e~ annihilations) (20) 

and we have not to worry about which variable was used. This will turn out to be important later 
on, where the experimental evidence of what has been conjectured so far is provided. 

Before considering this evidence, I would like to address what I consider a very important 
conceptual issue. Conservation laws at vertices are always expressed in terms of 4-vectorial quan- 
tities (like 4- momenta, tensors, etc). This is true even in the leading effect framework, where the 
conservation law reads: 
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On the contrary, when the experimental results are analysed, you always have to do with some 
scalar quantity like \fs or Mx (as the final product of any analysis usually is some scalar function 
of one variable) . 

The approach pursued here provides a consistent jump between the former and the latter as it 
allows to get from (21) the following scalar conservation law: 

Vs = A TO t + £tot (4 ') 

It seems to be a trivial step but, as previously noted, the world-wide use of Mx shows as it 
often happens that we do not recognize simple facts as a conservation law violation: the very 
use of M x in the framework of the effective energy violates the 4-momentum conservation law at 
interaction vertices and, consequently, of the whole interaction itself. 



6 Evidence of (n ch j vs A TO t universality 

In this section we show how it is possible to recover QCD universality by using Atot to plot data 
coming from different experiments. The proof is obtained by rescaling a result given in another 
paper [5]. 7 

It is appropriate here to briefly summarize the results of that work. It is a study of the 
quantity (n c h) in DIS processes. The main result is that (n c h)ms is more than 15% higher than 
(n c h) measured for p-p and e + e~ processes. This is well above the experimental uncertainties, so 
the conclusion was: DIS processes disagree with the leading effect phenomenology and bring to a 
QCD non-universality. Figure 2 shows the variable (n c h) as measured in [5]. 

The study was conducted in a reduced phase space. The charged tracks and the energies were 
only measured in the polar angle interval 9 LAB G [20°; 160°], see Figure 3. 

This cut is justified by a large improvement in the experimental resolutions and, as proved in 
the same paper by means of MC simulations, it does not affect the result in any way. The other 
cuts are not relevant and are those normally used to obtain an almost pure set of DIS events. 
More details may be found in [5] or below. 

Because of this cut in the polar angle, the variable Mx is substituted with M e g : it is again a 
hadronic invariant mass but this time it is not the total one; in fact, it is that measured in the A9 
considered. Again, MC simulations proved that the dependence between (n c h) and Mx (M e g) is 
not affected by the phase space reduction. 

From now on, unless otherwise stated, all variables refer to the reduced phase space without 
introducing any new notation to recall this fact. The only variable which notation will be changed 
is Mx to M e ff in order to help any comparison with [5]. 

Of course, the correction factor is the ratio between ^4tot and Mx and will be denoted by o. 
This factor will be applied to the abscissa in the plot showed in Figure 2, and will be obtained 
in a high energy approximation: this will cause a very slight mismatch at low energies even when 
Atot is used for re-plotting. This approximation could be avoided by re-analysing the data in 
terms of ^4tot- 

To give an estimation of the correction factor a that is suitable for our purposes, the general 
result (16) is useless as p™ D was not measured. Also that result is strictly valid only in the full 



Presently I have no access to row data, so I could not make a direct analysis. 
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Fig. 2 Figure taken from [5] showing the mean charged particle multiplicity in DIS events plotted as 
a function of M e s and confronted with the same quantity measured for other processes. The inner bars 
show the statistical errors, the outer ones show the statistical and systematic errors added in quadrature. 




Fig. 3 Typical DIS event in the ZEUS detector. The polar angle cut is shown: shaded areas are not 
considered in the analyses. Figure taken from [5]. 
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phase space. Thus the correction has been estimated by averaging over the two variables on which 
the said ratio depends, namely 9 and y (as showed below). The latter is a new variable that is 
defined in the next section and that expresses the unbalance in hadronic energy between the two 
hemispheres in a way that is suitable for our calculation. 

We proceed by estimatimating the quantity s ■ (A^ OT — M e 2 ff ) : 

s ■ [A% m M 2 S ] = [E^E 1 p"f] 2 [(E'f (p 1 ) 2 }^) 2 (P H ) 2 } = 

= (E H E 1 ) 2 + (p H p I cos9) 2 - 2E H E I p H p I cos 9- 

- (E H E') 2 - (ptp") 2 + (E'p") 2 + (£V) 2 = (22) 

= (pV) 2 (cos 2 6-1)- 2{E H E I p H p I ) cos 9+ 

+ {E I p H f + (E H p I f = ... 

where the index TOT has been suppressed and the absolute value of 3-vectors is denoted by 
suppressing the arrow above the vectors. 

Of course 9 is the angle between p 1 and p H that is, considering HERA kinematics, between 
p H and the z axis: 

9 = Z [p» ■ z) (23) 

This means that 9 is the polar angle, and will be identified with it from now on. Introducing the 
high energies approximation: 

p H ~E H 
p i ^ E i 

and substituting in (22) we get: 

. . . = (E H E I ) 2 (cos 2 9-1)- 2{E H E I f cos 9 + 2{E H E 1 ) 2 = 

= (E H E 1 ) 2 (1-2 cos 9 + cos 2 9) = (22') 
= [^^(l-costf)] 2 

To sum up: 

s ■ [A 2 TOT M 2 S ] ~ [E H E\1 - cos9)] 2 (E"' 1 ~ p^ 1 ) (24) 

In order to further simplify the last expression we are now forced to choose a reference frame. 
Again the best choice is the CMS, both to simplify s and for sake of visualization. Besides, it is in 
general the most "balanced" system, and this helps in avoiding as many troubles as possible with 
the high energies approximation. The only complication induced by this choice is that it forces to 
perform the needed transformations on data collected in the lab system. 

Calculating in the CMS we obtain: 

8 ■ [A 2 TOT - M 2 e ] ~ (2S*)L s (£ J )c MS (l - cos#) 2 CMS (25) 
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and using 

* = (C)c„ s (26) 
it becomes 

[^tot ~ M e 2 ff ] ~ (E H )l us (l - cos0) 2 CMS (27) 

As already mentioned, the ratio between ^tot and Mx is a function of two variables: AO and 
y. This means that, after having imposed the cuts used in the paper to be corrected, considering 
the collision parameter at HERA, and finally averaging over these two variables, we are left with 
a linear dependence of ^tot from M e g : 

Atoi oc M e g (fixed kinematical conditions) (28) 
. we have 

a(A6,y) (29) 



Thus, if in general we have 
Atot 



M eff 

the previous remark allows the calculation of the constant a — (a(A0, J/)), Ae y s that expresses the 
ratio between Atot and Mx at fixed kinematical conditions: 

a = (a(A6,y))^ Ae ^ (fixed kinematical conditions) (30) 
Note that a is invariant being the ratio between two invariants. 

The dependence of a from AO and y will be used later to show how it is possible to obtain 
a different value for a by tuning these parameters. This fact supports the correctness of the 
introduction of ^tot , as we obtain a value that allows DIS data to match with other processes 
curves if and only if the cuts used in [5] are taken into account. 

Using the notation introduced in (29), equation (27) can be rewritten as: 

[a 2 {A6,y)-l]-Mt s ~ (E H f cus (l - cos0) 2 CMS (27') 

or, isolating a, 



<r(AH. a) . L + ^P^-(l-cos^ MS (31) 



The relevant cuts and numerical values are: 

• #lab € [20°; 160°]: the already mentioned polar angle cut [5]. It must be re-evaluated in the 

CMS; 

• O^'Omin = ^ Gev, where the prime is used to refer to the final state and the "e" indicates 
the final state positron (this cut is imposed to guarantee an almost pure level of DIS events) . 
Even this value must be re-evaluated in the CMS; 

• /?boost = 0.935 is the value of the boost between the LAB and the CMS at HERA. 

and we are now ready to estimate the correction factor we are seeking, that is, the mean of a with 
respect to the variables it depends on: 



a = 



V <M.(r/(-E?)c„,)' 
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6.1 Average over the polar angle 

In this subsection only a prime is used to indicate CMS variables while unprimed varibles refer to 
the LAB system. Writing the appropriate transformations for the 4-momenta 

' E< = 7 (£7 - (3p z ) 

(33) 

p' z = j( Pz - f3E) 
and using again the high energies approximation in the form 
p z = p cos 9 ~ E cos 9 

(34) 

p' z =p'cos9' ~ E'cosO' 



we get 



E' = j(E-(3E cos 6>) 

(35) 

E'cosB' = j(E cos 9 - f3E) 



so that 



, ~/E{cos9 - f3) cos 9-/3 
C0S ° ~ 7 £(l-/3 cos 9) ~ l-/3cos0 (36) 

that is the well known formula for the aberration of light. The use of the last result is again cause 
of some overestimation of a at low energies. 

Working in the CMS allows to perform the mean over the polar angle by simply averaging over 
the distribution obtained in (31) as, if the sample is large enough, (0) cms ~ tt/2. 

Using the interval given in [5], 9 € [20°; 160°], and performing the needed transformations we get 
'20° — > 87.895° = 1.536 rad 

(37) 

160° — > 176.305° = 3.077 rad. 
These are the extremes of integration to be used. The integral yields 

/■3.077 

/ (l-cos0) 2 d6» 

((! - co ^) 2 )cms = 71 536 M - 2 -68 (38) 

so that 



/ 2 68 

a(y) = ( a(A9, y) ) (Ae) ^ ^1 + (Meff/£?g)2 (39) 

6.2 Average over the new unbalance variable "y" 

As anticipated above, we now define a variable that suitably represents the unbalance in hadronic 
energy between the two hemispheres. We call this variable y and define it as follows: 

(40) 

E» e = {l-y)-E» 
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that is, y is the fraction of hadronic energy measured in hemisphere 1 (that of the "target region", 
where the proton remnant is detected). 

Using y we can express M x as a function of E"£% in a very useful way. In fact, assuming that 
~ — k^ (where k H represents the hadronic 3- momentum versor) we have 

M x ee v/feF = [(gf f + (gf f + 2E?E? 2p» ■ p ?] 1/2 ~ 

(41) 



that together with definition (40) yields 



M X c <l^E?E* = 2EZ°Vy~(l=y) (42) 

We remark that this rescaling method is applicable only if we have no more than 2 hadronic 
jets. In case of 3 or more jets, another method should be used or a direct data analysis must be 
performed. 

The hadronic energies we are considering are "real" hadronic energies, that is, the energies 
effectively produced at the interaction vertex, regardless of the cuts applied to perform the analysis. 

This does not create problems in case of a phase space reduction as there is no dependence of 
Mx (-E"o°) from the phase space [5]. On the contrary, the cut in the final positron energy must 
be considered, as the events that do not meet this condition do not become part of the statistics. 

As (y) CMS = 1/2 holds, for the same reason exposed in the previous section, it is sufficient to 
calculate the mean of the distribution obtained in the (42). Had there been no cuts in the leading 
particles energy, this integral would have yielded 7r/8, that is half the area of a circle of radius 
1/2. This is not our case and we will have to take into account the relevant cuts, one of which 
has been already mentioned and is: 



(20£ N B = 8Gev 
that, as usual, must be re-evaluated in the CMS: 



I \ LAB 




(43) 



= (K)™-\I^4 * 5.467 (E' e )™ ~ 43.7 Gev 



where we have: 

• considered that e z are antiparallel, from where the change in sign at the second equality 

• used the high energy approximation and the strong collinearity of the e^ UT (again at the 
second equality) 

• used the value of /5 BOO st for /? 

The highest leading longitudinal momentum in the target region is detected when a proton 
takes the role of the leading particle: as this maximum value, we use that representing the limit 
between the leading physics and the diffractive physics (in the sense described in [1] or [6]), namely: 

TO Wx = (x F )™ ■ = 0.8 £1- (44) 
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where xf is the Feynman variable that represents, in practice, the fractional longitudinal momen- 
tum of the leading particle. 

It should be noted that the relation between energies or momenta expressed in fractional terms 
are approximately invariant. In fact: 

'E' hEAD =~f(E L -p.p L ) 



(45) 



E' L = j(E L -(3-p L ) _ 7 (1-/?) & 
E'j rtE'-tf-f 1 ) 7(1-/3)'^ 



E 1 



so that the energy of a particle expressed with the Feynman variable is about the same in two 
different frames at high energies and in a collinear approximation. The leading particle, as such, 
respects these approximations. 

The two values just obtained allow to find the largest interval in hadronic energy compatible 
with the cuts: 



( F H\cms = 4 V* 

V I /MIN o CO 



^2 Jmax — o V^e^MiN 



30.0 Gev 



106.5 Gev 



from which the minimum for y immediately follows: 



V:, 



/CMS 



CMS 



0.220 



(46) 



In order to give an estimation for y MAX , it is necessary to make some assumptions regarding 
x F m and • A good choice, in view of a comparison with ISR data, is to adopt the same 

value of x F 1N chosen at ISR, namely 0.4 : 



( ry, \CMS r. A TJiINC 

It will also be assumed that the final state positron loses at least the 30% of its energy, value 
above which the diffractive events begin to prevail over the DIS ones: 



(K) 



I ^CMS 
MAX 



0.7- 



105 Gev 



From the previous two equations we get 
3 



(En 



H \ MAX 

CMS 



H\MIN 

CMS 



(E?) 
that in turn yields 



5 2 



90.1 Gev 
(E' e )Zl * 45.0 Gev 



(47) 



(En 



H\MAX 

CMS 



(En^+(En 



MIN 

CMS 



~ 0.667 



(48) 
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Summarizing: 



(En 



CMS 



Umax — 



(En 



MAX 
CMS 



pH \ MAX 
/CMS 



0.220 



~ 0.667 



(49) 



(£f )«- + (E?)™Z 

and it is now possible to give a numerical estimation for the mean of the distribution in hadronic 
energy a depends on: 

pO.667 



1ERA 
' CMS 



/•U.DD I 

/ Vv( 

JO. 220 



0.667- 0.220 



or, in our case, 

(m x )z r ; = 2Ko T d -(v^t^)):: 



1 - y) dy 



~ 0.48 (cuts considered) 



0.96 • E" 



(50) 
(51) 



7 Final result 

Given the previous results and considerations, we obtain the following correction factor a to be 
applied in abscissa in [5]: 



a = 



((1- cos 0)2 ) / 2.68 , « 

1.98 (52) 



\ (M eff /(^) CMS ) 2 V (°' 96 ) 



So, the plot will be showing the (n c h) vs Atot dependence if we multiply DIS data abscissas by 
a factor a= 1.98. We recall once more that e + e~ data do not necessitate any rescaling while for 
p-p data an estimation of the correction factor is given below, but we can anticipate that it is 
negligible. The result is showed in Figure 4. 

It is evident that, particularly at high energies, DIS data and other processes data follow the 
same curve. The quantity (n c h) is again a Universality Feature. 



8 p-p data rescaling 

Recalling that the ISR was balanced and that the full phase space was considered, we can conclude 
that the analyses performed using E"*° (the almost totality) were automatically made in terms 
of Atot (see equation 11). This means that those works do not necessitate any correction. 

On the contrary, in some other works, the variable Mx was employed [6, 13, 14, 17]: these 
do need to be corrected for the reasons exposed in the introductory section. This may seem in 
contrast with the fact that no differences were found at that time when data were analysed both 
with -E"ot and Mx but it is not so, the reason being in the low precisions and energies then 
achievable. 

As an example, we estimate the difference between Atot and Mx that could have been revealed 
in the early works at ISR where the only condition was the presence of a leading proton leaving 
the vertex with 



x F e [0.4 , 0.8] 



(53) 



adp header will be provided by the publisher 



19 



18 



v 



-|g _ ■ ZEUS 95 Prelim., inclusive 

▼ ZEUS 94-97, current region 
(multiplied by a factor of 2) 



14 



12 



10 - 



8 - 



• Rescaled ZEUS 95 data 




4 - A A A 



2 - 



A_ 

Ia a 



b eVTASSO 
o e + e" PLUTO 
■k eVJADE 
a e + e MARK I 
• pp ISR 



10 



'TOT 



(GeV) 



10 



Fig. 4 The same figure as 2 but with DIS data (ZEUS 95) rescaled by the factor a ~ 2 and plotted in 
violet: now the X-axis variable is ^4tot- The universality of (n c h) is evident, particularly at high energies 
(a being an asymptotic value). 



from which it follows that 



MAX — ^ MIN 



(54) 



and this in turn allows to get the maximum value of a in the "worst" situation: 

^TOT (-EtotOcms.lab 



Mv 



My. 



("worst" situation) 
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where again the high energies approximation has been used in the second step. Thus, in a single 
event and in the worst situation, the difference between ^4tot and Mx could had been as large 
as 15%, but this is a limiting event; when averaging over all the events the difference becomes <C 
15% and is consequently hidden by the large uncertainties that affected {n c h) at those times. It is 
therefore possible to safely employ ISR p-p data analysed in terms of M x to make a comparison 
with other processes data if a very high precision is not required. 



9 Dependence of a on A.0 and A.y 

In this last section the dependence of a on A9 and Ay is discussed with the main purpose of 
showing how it is improbable for the result obtained in this paper to be fortuitous. 

In fact, changing kinematical conditions (using other values for A9 and for Ay) would bring to 
different values of a and the agreement between e + e~, p-p and DIS achieved in this paper could 
have not been obtained without taking a full account of the correct kinematical conditions. 

Equation (31) shows that a increases with 8 and decreases with the width of the interval in y 
(centered at 0.5). We examine the situations when a is as high or as low as possible. We also 
check what happens in case of no phase space reduction. 

It is of course safer to vary AO than Ay as from the latter depends the purity of the DIS 
sample. Only safe values for the minimum and the maximum of y will be considered. 

A maximum for 2(y / y(l — y) ) is obtained by selecting a sample of balanced events: 

y=\ — 2(V^T^)) max = 1 (56) 

while unbalanced events induce a minimum for a. If events with a fast leading proton and a low 
energy outgoing positron are chosen (events unbalanced in E H towards hemisphere 2) 

' x F € [0.6; 0.8] 

(57) 

JWn S € [43,7; 70] Gev 
then, in the same way used earlier, we get: 
'y MIN = 0.220 

(58) 

= 0.431 
from which 

(Vv0^y)) Mm = 0-464 (59) 

a value that may be used to re-evaluate a under modified conditions. As showed in (18') the 
analytical minimum for a is 1: 

fl =Me7^ (60) 

but it is not possible to reach this value in true experiments as measurements should be performed 
in very small and collinear intervals in 9 (say 9 LAB <G [0°; 5°]). Perhaps, a realistic possibility could 



adp header will be provided by the publisher 



21 



be 8 CMS £ [0°; 90°] that is 9 LAB £ [0°; 20°] (a choice that implies a poor resolution anyway). In such 
a case 

<(l-costf) 2 CMS )^0.23 (61) 
which yields 



23 

a MIN ~ Jl + — ~ 1.1 (62) 

that is, in these conditions Atot and Mx would be undistinguishable in practice, again because of 
poor resolutions. Such a study would have never been able to reveal any difference between ^4tot 
and M e g and the universality features would have been evident even using M e g as it happened in 
some previous experiments. 

In order to get much greater values for a than that obtained in this work, we should go to the 
other side of the polar angle range, imposing, for example, ^cms G [150°; 180°]. But this would 
mean to say 9 LAB £ [174°; 180°], again a ridiculous interval. We are thus forced to conclude that 

a MAX >2 (using (y/y(l - y) ) m _ n = 0.464) (63) 

and note that this study has already been performed in in a kinematical region where a reachs or 
is very near to its maximum. 

Finally, it is interesting to see what would have been obtained by working in the full phase 
space. Imposing 

A0 LAB = A0 CUS = [0°; 180°] (64) 

we get 

((l-cos0) 2 CMS ) = 1.5 (65) 

so that 



/ 1.5 

a (FULL PHASE SPACE) W 1 ~l~ (Q Q6)2 ~ 1-62 (66) 

a factor implying a difference that is about 40% smaller than what has been found in the reduced 
phase space. This result, combined with the worse resolution at the two ends of the ZEUS 
calorimeter, shows that a full phase space study would have been loosely able to reveal that 
smaller difference between A^ot and Mx . 

This discussion shows how the value estimated for a and the consequent gathering toward the 
same curve of data points related to different processes, has a quite pronounced dependence from 
the "contour" conditions. As an example, using a( FULL PHA se space) to rescale DIS points in Fig. 4, 
would keep the mismatch between DIS and other data points evident. 

This implies that a possible fortuitousness of the agreement obtained can be regarded as quite 
unlikely. It is of course necessary to perform direct analyses to check if Atot works well in 
different conditions and kinds of experiments, particularly in those involving the indirect leading 
effect: these processes in particular appear to be the best candidates to check the geometrical 
justification for introducing A^qt ■ 
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10 Conclusions 

After reconsidering under the correct geometrical viewpoint the kinematics of the leading effect, 
a new invariant has been introduced. It is likely to be the universal Lorentz-invariant quantity 
effectively available to multiparticle production and, as such, the variable that correctly describes 
the hadronization processes in any kind of particle collisions. Consequently, it has been called the 
Total available quadri-scalar and denoted with Atot ■ 

The universality of (n c h) has thus been recovered after that other works [5] had claimed its 
non-universality when comparing DIS with other kinds of processes and analysing in terms of the 
invariant hadronic mass Mx- 

This result definitely rules out Mx as the quantity representing the so-called "hadronic energy" 
and supports Atot as its most reasonable successor. In order to exclude other yet-to-be-proposed 
possibilities to solve the universality problem, and eventually adopt ^4tot to study and describe 
available and future experimental results, further analyses concerning other physical quantities 
need to be performed. 

Appendix A: 

A comparison with the generalized parton distribution approach 

Finally, I would like to make a comparison and discuss similarities and differences between this 
work and the generalized parton distribution (hereafter: GPD) approach to the physics of particles 
collisions. For a recent review on the topic see [22] and references therein. I have found some 
similarities in the two viewpoints, but the main difference stays: the incorrect use of invariant 
masses to study particles collisions in general and hadronization in particular. 

Of course the GPD formalism digs deep into parton dynamics and enables to perform explicit 
calculations concerning exclusive processes and to study some (otherwise unaccessible) quantities 
like, for instance, 3-dimensional (transverse) distributions of partons inside the target hadron (see 
Impact Parameter Representation, Sect. 3.10 of [22]) or helicities. It allows to derive interesting 
information about the internal structure of hadrons as well as to calculate form factors for sub- 
processes whose experimental study is very difficult or even not feasible, like graviton-parton 
interactions or gluonic currents. 

The present work is based on 4-dimensional geometrical-kinematical aspects of the collision 
processes and consequently distinguishes longitudinal and transverse components of involved 4- 
momenta, just as the GPD formalism does. This surely is a common point to both approaches 
but general results and, in particular, aims are quite different. 

In fact the present work concerns with a different aspect of particle reactions, namely the global 
and inclusive properties of multi-hadronic final states, which must show a global universality as a 
consequence of QCD universality, as has been extensively discussed in the paper. 

It does not discuss exclusive processes and does not allow to compute associated amplitudes 
(even though it tells which is the correct variable to use to perform such a task) , but concentrates 
instead on inclusive features of various types of reactions. We describe the collision as a 2-step 
process, and this is another similarity with GPD approach, but there is a main difference: our 2 
steps are "visible" in the final states: the leading particle (or target remnant) is not considered 
as a part of the multi-hadronic final state. After this subtraction, any feature one could think of 
should be the same, independently from reacting particles or nominal energies. I have showed this 
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is indeed the case for the most essential quantity related to the study and description of global 
properties of final states, i.e. (n c h) (the mean charged particle multiplicity). 

The GPD approach is unable to recover this universality, the reason being the use of invariant 
masses. It forces to account for different kinematics when comparing processes of different nature 
(see e.g. Sect. 6.5.1 of [22]), when instead there is no difference at all. There are no "diffferent 
kinematics", we are simply looking at things from the wrong point of view (i.e. using the wrong 
variables) . 

The use made of invariant-mass-type quantities, like Mx (which does not take into account the 
transverse contribution from the leading particle) or W 2 (which, even worse, does not take into 
account the leading particle at all) will never allow to correctly analyse any kind of experimental 
data. 

Of course GPD could provide very detailed predictions, were the distinction between the leading 
particle and the rest of the final state implemented into the formalism. It would be very interesting 
to use the power of GPD to extrapolate results to areas not accessible to the experiment in order 
to check and extend the study of hadronization under the ^4tot framework toward those events 
where the leading particle is not detectable (and thus its subtraction becomes difficult if not 
impossible). 

At any rate, it is of fundamental importance to abandone at once - both for the reason presented 
here and for experimental evidences of non-universality (also presented and discussed in the paper) 
- the use of variables like s , Mx , W 2 and similar when their use is evidently improper. 

I suspect that many controversial results obtained in the small-x regime (where the leading effect 
mostly plays its role: high longitudinal momentum of the leading particle: higher subtraction), 
as, for instance, those cited in Section 8 of [22], could be settled were Atot and related variables 
used to re-analyse available or new data. 

If the use of Atot (and its derived or derivable quantities) is to be implemented in the QCD 
description of both hard partonic subprocesses and the hadronization phase, I believe that the 
longitudinal-transverse separation-oriented GPD formalism, showing the discussed similarities 
with the 'Atot approach", could facilitate efforts toward this scope. 
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